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We present a theoretical approach to calculate the molecular magnetic anisotropy parameters, Dm 
and Em for single molecule magnets in any eigenstate of the exchange Hamiltonian, treating the 
anisotropy Hamiltonian as a perturbation. Neglecting inter-site dipolar interactions, we calculate 
molecular magnetic anisotropy in a given total spin state from the known single-ion anisotropics of 
the transition metal centers. The method is applied to Mn^Ac and Feg in their ground and first 
few excited eigenstates, as an illustration. We have also studied the effect of orientation of local 
anisotropics on the molecular anisotropy in various eigenstates of the exchange Hamiltonian. We 
find that, in case of Mn^Ac, the molecular anisotropy depends strongly on the orientation of the 
local anisotropies and the spin of the state. The Dm value of Mn\iAc is almost independent of the 
orientation of the local anisotropy of the core Mn(IV) ions. In the case of Fes, the dependence of 
molecular anisotropy on the spin of the state in question is weaker. 

PACS numbers: 75.50.Xx, 75.30.Gw 



I. INTRODUCTION 

Following the synthesis and the discovery of exotic 
properties such as quantum resonant tunneling (QRT) 
in the single molecule magnet (SMM) Mn^Ac during 
the 1990s, there has been a flurry of activity in the field 
of molecular magnetism [TJ EJ |3j |4j . This has led to the 
synthesis of new systems such as Feg as well as to the 
observation of new phenomena such as quantum coher- 
ence SMMs are mainly high nuclearity transition 
metal complexes with a high-spin ground state (Sgs)- 
They are also characterized by large uniaxial magnetic 
anisotropy [TJ. The Hamiltonian corresponding to the 
magnetic anisotropy of a molecular system can be writ- 
ten as, 

H D = S M ■ V {M) ■ S M (1) 

where, Sm is the spin operator for the total spin of the 
molecule and I?( M ) is the magnetic anisotropy tensor of 
the molecule. In usual practice, the anisotropy tensor is 
diagonalized and the principle axis of the molecule would 
correspond to the eigenvectors of the tensor. Since in 
most physical situations, the quantity of interest is the 
energy gaps between the otherwise degenerate states split 
by the magnetic anisotropy, the condition of zero trace 
is imposed on the X>( M ) tensor. If, D™ x , D™ Y and Df z 
are the molecular anisotropies along the three principal 
directions such that D xx + Dyy + D zz = 0, we can 
define two parameters, Dm and Em given by, 

Dm = Df z - \ (Dxx + D¥r) 

E m = - (Dxx - Dyy) 
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where, Dm and Em are called the axial and rhombic 
anisotropies respectively. This leads to the common form 
of the magnetic anisotropy Hamiltonian of a SMM, 

H M = D M [s 2 z - ^S(S + 1)^ + E M (s 2 x - &) (2) 

For the single molecule magnet to have nonzero mag- 
netization in the ground state, it is necessary that the 
anisotropy constant Dm in the spin Hamiltonian (Eq. [2]) 
of the complex be negative; this ensures that the ground 
state of the system then would correspond to the high- 
est magnetization state of the molecule, in its high-spin 
ground state. This requirement of negative Dm, besides 
a high-spin ground state makes it hard to tailor SMMs. 
The second order transverse or rhombic anisotropy given 
by the last term in Eq. [2] allows transition between states 
with spin S that differ in their M s values by two. Em will 
be zero if the S x — S Y operator does not remain invari- 
ant under symmetry of the molecule. In this case higher 
order spin-spin interaction terms are required to observe 
QRT. For example, the D 2 d symmetry in Mn^Ac pro- 
hibits the existence of first order rhombic anisotropy. 
Thus, the parameters Dm and Em govern the quantum 
tunneling properties of a SMM and inputs from theoreti- 
cal modeling could help in designing the architecture for 
the synthesis of SMMs. 

Theoretical modeling of SMMs presents two difficul- 
ties. Firstly, the complexes contain many spin centers, 
and often these centers have different spins as in the case 
of Mn^Ac, where the four Mn(IV) ions have spin-3/2, 
while the eight Mn(III) ions have spin- 2. In these sys- 
tems, usually there exist multiple exchange pathways be- 
tween any given pair of ions leading to uncertain mag- 
nitude and sign of the magnetic interactions in the sys- 
tem. The topology of magnetic interactions also is often 
such as to result in magnetic frustration leading to closely 
spaced low-lying states of unpredictable total spin. Thus, 



FIG. 1: Schematic of possible exchange interactions in 
Mn\iAc SMM. The peripheral Mn(III) ions represented by 
blue circles correspond to spin-2 sites and those represented 
by yellow circles are the core Mn(IV) ions each of spin-3/2. 
Js are the strength of superexchange interaction with Ji=215 
K, J 2 =J 3 =85.6 K, J 4 =-64.5 K 0. 



solving even the simple Heisenberg exchange Hamilto- 
nian of these systems turns out to be a challenge. If we 
do not employ a spin adapted basis to set-up the Hamil- 
tonian matrix we could encounter convergence difficulties 
associated with closely spaced eigenvalues even when the 
corresponding cigenstatcs belong to different total spin 
sectors. However, the assorted spin cluster that a SMM is 
renders construction of spin adapted basis difficult. This 
problem has been addressed by resorting to a valence 
bond scheme for construction of the spin adapted basis 
and we are now in a position to block-diagonalize the 
Hamiltonian by exploiting both spin and spatial symme- 
tries [8]. Thus, while the challenge of accurately solving 
the exchange Hamiltonian of a large assorted spin sys- 
tem with arbitrary topology of exchange interactions is 
within grasp, the challenge of computing the magnetic 
anisotropy constants of a SMM still remains. 

The magnetic anisotropy in an isolated ion arises from 
explicit dipolar interactions between the unpaired elec- 
trons in the magnetic center as well as from relativistic 
(spin-orbit) interactions. The former is the main ori- 
gin of zero field splitting observed in triplet states of 
conjugated organic molecules such as naphthalene [10] , 
However, in systems containing heavier elements the rel- 
ativistic effects dominate. In a system with several mag- 
netic centers such as a SMM, given the spin and sin- 
gle ion anisotropy of each magnetic center, the mag- 
netic anisotropy could arise both from dipolar interac- 
tions between magnetic centers and relativistic effects. 
Theoretically, magnetic anisotropy of an isolated mag- 
netic center in an appropriate ligand environment can 
be computed within density functional techniques (DFT) 
[TT| [12]. However, in the case of the SMMs with many 
magnetic centers, such a calculation is both conceptu- 
ally and computationally difficult, since DFT methods 



FIG. 2: Schematic of local (a;, y, z) and laboratory (X, Y, Z) 
coordinate axes in MnviAc. The blue, green and red spheres 
correspond to Mn(III) (spin-2), Mn(IV) (spin-3/2) and oxy- 
gen ions respectively. The arrows indicate the Mn-O bonds 
on which the chosen local x-axis has maximum projection. 

do not conserve total spin. The usual approach in these 
cases is to carry out simple tensoral summation of the 
anisotropies of the constituent magnetic centers to ob- 
tain the magnetic anisotropy in the SMMs along the lines 
of an oriented gas model employed in the calculation of 
macroscopic nonlinear optic (NLO) coefficients from iso- 
lated molecular NLO coefficients [121 HH US]. Such an 
approach, in the case of SMMs, suffers from the draw- 
back that the anisotropy constants so computed are in- 
dependent of the total spin state of the molecule. In 
this paper, we treat the anisotropic magnetic interaction 
between the magnetic centers as a perturbation over the 
magnetic exchange Hamiltonian of the SMM. We exactly 
solve for the various total spin states of the unperturbed 
Hamiltonian and from the desired eigenstates obtain the 
spin-spin correlation functions necessary to compute the 
anisotropy constants. In the next section we describe the 
method in detail. In the third section we present the re- 
sults of our studies on the two SMMs, Mn^Ac and Fes- 
In the fourth section we summarize our studies and out- 
line the extension of this method for the computation of 
higher order magnetic anisotropy constants. 

II. FORMULATION OF THE METHOD 

We treat the exchange Hamiltonian between magnetic 
centers in the SMMs as the unperturbed Hamiltonian, 

Hq = 2^ Jij&i ■ Sj (3) 

(ij) 

where, (ij) runs over all pairs of centers in the model 
for which the exchange constant is nonzero, Si is the spin 
on the ith magnetic center. In SMMs such as Mny^Ac 
the spins at all the magnetic centers are not the same and 





FIG. 3: Schematic diagram showing the directions of local 
anisotropy in Mn\iAc. The single-ion anisotropics of all the 
Mn ions are directed along the laboratory Z axis (Scheme 1). 



the exchange interactions are shown in Fig. [T] H can 
be solved exactly for a few low-lying states in a chosen 
spin sector by using methods that have been described 
in detail elsewhere [9]. 

The general anisotropic interactions in a collection of 
magnetic centers is treated as a perturbation with Hamil- 
tonian H[ is given by, 



^ = ~££££^ 



ij,af3 



(4) 







where, the indices i and j run over all the magnetic cen- 
ters and a and (3 run over x, y and z directions of the 
ion. The contributions to inter-center anisotropy con- 
stant arise due to dipolar interaction between the spins 
on the two centers as well as due to relativistic effects. 
In the former, Dij^p is given by, 
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ij,a/3 
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(5) 



where IZij (Rij) is the vector (distance) between the mag- 
netic centers i and j, g is the gyromagnetic ratio and [Ib 
is the electronic Bohr magneton; the expectation value 
in Eq. [5]is obtained by integration over spatial coor- 
dinates ]l6j . Approximating the expectation values of 
the distances by the equilibrium distances, the Dij a p in 
Eq. [5] and by computing the necessary spin-spin corre- 
lation functions, we can obtain the molecular v£p ten- 
sor |17j . The eigenvalues of this matrix give the princi- 
pal anisotropy values and imposing the condition of zero 
trace of the matrix yields molecular magnetic anisotropy 
constants due to spin-spin interactions. Our computation 
of the magnetic anisotropy constants, assuming only spin 
dipolar interactions for the SMMs MnnAc and Fes gives 
negligible values of the anisotropy constants compared to 
the experimental values of Dm = -0.7 K and -0.28 K re- 
spectively in the S=10 ground state (TSJ HH]- Hence in 
what follows, we completely neglect the contribution of 
spin-dipolar interactions and focus only on the magnetic 
anisotropy of the SMMs arising from the anisotropy of in- 
dividual magnetic centers. The latter is a consequence of 



FIG. 4: Schematic diagram showing the directions of local 
anisotropy in Mn\iAc. The z-component of the single-ion 
anisotropics of all the Mn(III) ions are inclined at an angle 
6 to the laboratory Z and while that of the Mn(IV) ions are 
kept fixed at ~ 48° (Scheme 2). 



mainly spin-orbit interactions. We now assume that the 
interactions responsible for magnetic anisotropy are short 
ranged and neglect inter-center contributions to magnetic 
anisotropy in Eq. [4] This is justified since relativis- 
tic (spin-orbit) interactions, largely responsible for the 
anisotropy is short ranged (falling off as l/r 3 ) and the 
distances between the magnetic centers is much larger 
compared to the ionic radius of the transition metal ion. 
The resulting perturbation term is given by, 



Hi = ££ £ Di^pS?^ 



(6) 



where, only on-site terms are retained. If the individual 
magnetic centers have different principal axes then we 
choose a laboratory frame and project the local tensor 
components on to the laboratory frame. In such a case, 
Eq. [6] gets modified to, 

Hi = £ £ £ £ £ CuaC^mpDi^pSfS^ (7) 



Where, Cn a are the direction cosines of the local axis 
of the ith magnetic center with the l(m) being the co- 
ordinate of the laboratory frame and a((3) being the lo- 
cal coordinates. Since, the Hamiltonians in Eq. [I] and 
[7] are equivalent, we can equate the matrix elements 
(n,S Jf ,M|Ai|n,SA f ,M') and (S M ,M\H M \S M ,M') for 
any pair of eigenstates of the exchange Hamiltonian in 
Eq. [3] I M) and \M') correspond to a state n with spin 
Sm in which we are interested. Calculating these ma- 
trix elements for Hm is straightforward from the algebra 
of spin operators. However, evaluation of these matrix 
elements between eigenstates of Hq, requires computing 
them in the basis of the spin orientation of the sites. From 
a given eigenstate of Ho, \n,S,M), we can compute all 
eigenstates with different M values by using ladder oper- 
ators corresponding to spin S. Given a S value, the above 
condition would give rise to (2S + l) 2 equations, while 
the tensor T>^ M \ has only nine components. Thus, for 
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FIG. 5: Schematic diagram showing the directions of local 
anisotropy in Mn^Ac. The z-component of the single-ion 
anisotropies of all the Mn ions are directed along the plane 
perpendicular to the laboratory Z axis (Scheme 3). 



the Mn 12 Ac system, with ground state spin of 10, there 
would be 441 equations and we have more equations than 
unknowns. However, we could take any nine equations 
and solve for the components of the tensor 2?( M ) and we 
would get unique values of the components. This is guar- 
anteed by the Wigner-Eckart theorem and we have also 
verified this by solving for the 2?( M ) tensor from several 
arbitrarily different selections of the nine equations. 



III. RESULTS AND DISCUSSION 

We have computed Dm and Em values for both 
Mn\iAc and Feg systems for different orientations of 
local anisotropy. We have used the single-ion anisotropy 
values quoted in the literature for complexes of these ions 
in similar ligand environments [IB ED HH US- While 
discussing the results, we refer to the local axis of the 
ions as x, y and z and the laboratory axis is denoted as 
X, Y and Z. The laboratory frame we choose can be 
arbitrary. This is because, on determining in the 

laboratory frame, we diagonalize it and the principal axis 
of the molecule is given by the eigenvectors of the X>( M ) 
matrix. The principal axes of the molecule are unique 
and do not depend on the laboratory frame that is se- 
lected. We have computed the anisotropy parameters for 
both these systems as a function of the angle 9 which the 
z-axis of the ion makes with the laboratory Z-axis. The 
orientation of z-component of the single-ion anisotropy 
in every site is shown in Fig. 
and 3) for Mn^Ac and in Fig. 
and 6) for Feg, systems respectively. Once the z-axis (z) 
of the ion is fixed, then x is obtained by Gram-Schmidt 
orthogonalization procedure. Though the choice of this 
vector is arbitrary in a plane perpendicular to z-axis, we 
have fixed the direction of x such as to have maximum 
projection along a M-0 (M=Mn, Fe) bond in Mn^Ac 
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FIG. 6: Variation of Dm as a function of 8, the angle the 
z-component of local anisotropy of Mn(III) ions makes with 
the laboratory Z-axis for scheme 2, in eigenstates with total 
spin 10, 9 and 8. The orientation of Mn(IV) ions is kept 
fixed at ~ 48° from the molecular Z-axis. The curve with 
filled circles correspond to the variation of Dm, when the 
local anisotropies of the core Mn(IV) ions only are rotated 
and those of Mn(III) ions are fixed along the Z axis. 



10 1. If O is the vector 



as well as in Fe§ (Fig. [ijand Fig. 
connecting a M site and a neighbouring O ion, then we 
obtain x from, 



0-0 



(8) 



Then, the y-axis of the ion is obtained by taking cross 
product of the z and x, y — z x x. These three mutually 
orthogonal vectors are then normalized to obtain the or- 
thonormal set of coordinate axes x, y and z of the ion 
centre. The single ion local axes is represented in the 
laboratory frame as, 



X = 


Ci,xxX - 


- Ci^YxY - 


1- Ci^zxZ 


y = 


Ci,XyX - 






z = 


Ci : xzX - 


- Ci, Yz Y H 





(9) 

where, Cs are the direction cosines of Eq. [7] and the in- 
dex i correspond to the site i. The procedure is repeated 
for every magnetic ion to obtain coordinate axes set x, y 
and z and the direction cosines in each case. We have ob- 
tained the magnetic anisotropy parameters Dm and Em 
for Mni2^4c and Fes clusters as a function of the angle 
of rotation of the local z-axis with respect to the labora- 
tory Z-axis. In the following subsections, we discuss the 
results for the two clusters. 



A. Magnetic anisotropy in Mn^Ac SMM 

We have first obtained the ground state and few ex- 
cited states of the Mn\i Ac system by exactly solving the 
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TABLE I: Dm values of ground and excited states of Mn^Ac 
under various schemes in K. For scheme 2, we have presented 
the Dm values only for 0=26.2° for which the Dm value of 
the ground state matches with the experimentally observed 
value. 



State 


D M (K) 


Scheme 1 


Scheme 2 


Scheme 3 


Ground state (S=10) 


-0.8028 


-0.7209 


0.3991 


First excited state 
(S=9) £ s =35.1 K 


-0.6722 


-0.6105 


0.3341 


Second excited state 
(S=8) £ 9 =60.4 K 


-0.5009 


-0.4664 


0.2488 



unperturbed Hamiltonian given in Eq. [3] using the ex- 
change interactions shown in Fig. [I] [3] . The ground state 
of the system corresponds to total spin 10 with a total 
spin 9 excited state at 35.1 K from the ground state. The 
second excited state occurs at 60.4 K from the ground 
state and corresponds to total spin 8. To obtain the 
molecular anisotropy values in these cigcnstatcs, we have 
used different single-ion axial anisotropy values of -5.35 
K and 1.226 K respectively for Mn(III) and Mn(IV) 
ions. We have also introduced transverse anisotropy of 
0.022 K and 0.043 K for Mn(III) and Mn(IV) sites 
respectively. We have studied the variation of molecu- 
lar anisotropy as a function of orientation of the local 
anisotropies by rotating the local T> tensor around the 
molecular Z-axis. Scheme 1 shown in Fig. [3] corresponds 
to the case when all the single-ion z axes are pointed par- 
allel to the laboratory Z direction. In scheme 2 (Fig. |5|, 
we have fixed the orientation of the local anisotropics of 
the core Mn(IV) ions along the line joining the ion and 
the molecular centre (~ 48° from the laboratory Z-axis) , 
while the anisotropies of the Mn(III) ions is rotated 
and the angle which it makes with the molecular Z-axis 
is defined as 9 (refer Fig. [4]). The orientation of the local 
anisotropy of Mn(III) ions for which we get the best 
agreement with experiments corresponds to 9 — 26°. In 
scheme 3 (Fig. [5| , we have restricted the z-component of 
single-ion anisotropy to the laboratory X — Y plane. We 
have studied the variation in the molecular anisotropies 
in these schemes for ground and the excited eigenstates. 
We show in Table [TJ the Dm values for the ground and the 
excited spin states of the molecule for schemes 1, 2 and 3. 
We note that when the local anisotropies are systemati- 
cally varied, there is a very large variation in the molecu- 
lar anisotropy as a function of the local orientation (Fig. 
[6]). This seems to be true for all the states of Mn^Ac 
we have studied. We note that given the orientations of 
local anisotropies, the actual molecular anisotropy values 
are different in different spin eigenstates. This may be 
rationalized from the fact that the energy gaps between 
the ground and the excited states are large as a conse- 
quence of which the spin correlations in these states are 
very different. We also note that in all cases, from the 
eigenvectors of the 2?( M ) matrix, we find that the choice 
of our laboratory frame is very close to the principal axis 




FIG. 7: Schematic of exchange interactions in Fes SMM. Js 
are the strength of superexchange interaction with Ji = 150 K, 
J 2 =25 K, J 3 =30 K, J 4 =50 K [9]. 



z 




FIG. 8: Schematic diagram showing the directions of lo- 
cal anisotropy in Feg. The single-ion anisotropies of all the 
Fe(III) ions are directed along the laboratory Z axis (Scheme 
4). 



of the molecular system. 

We also examined the role of magnetic orientations of 
the core Mn(IV) ions (s=3/2) and the crown Mn(III) 
ions (s=2) in determining the molecular anisotropies by 
fixing the single ion orientation of the crown Mn(III) 
ions at 0° and rotating only the orientation of the core 
M n(IV) ions systematically. The variation of Dm for the 
S=10 ground state as a function of rotation of the local 
anisotropies of the core Mn(IV) ions is shown in Fig. [6] 
We find that the molecular anisotropy is not sensitive to 
the local orientations of the core Mn(IV) ions while the 
orientation of the crown Mn(III) ions control the vari- 
ation of the molecular magnetic anisotropy in Mn^Ac. 
It should be noted that, in the case of Mn^Ac, Em van- 
ishes by virtue of the D 2 d point group symmetry to which 
the molecule belongs. 
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B. Magnetic anisotropy in Fes SMM 



We have also computed the values of molecular 
anisotropy for the -Fes molecular magnet. The unper- 
turbed Hamiltonian in Eq. [3] is exactly solved using 
exchange parameters, Ji=150 K, J 2 —25 K, J 3 =30 K, 
Ji=50 K (Fig. [7]) [3]. The ground state of the system cor- 
responds to total spin S=10 with a S=9 state at 13.56 K, 
a S=9 state at 27.28 K and a S=8 state at 28.33 K above 
the ground state. To calculate the magnetic anisotropy 
of Fes, w ^ have taken the single ion axial and rhombic 
anisotropy values for Fe(III) centers to be 1.96 K and 
0.008 K respectively. Using these, we have computed the 
molecular anisotropy values for three schemes (Schemes 
4, 5 and 6) shown in Fig. [8] [9] and 10 Scheme 4 corre- 
sponds to the case wherein the single-ion anisotropy of all 
the Fe(III) ions are pointed along the laboratory Z di- 
rection. In scheme 5, the anisotropies of the Fe(III) ions 
are inclined at an angle 9 to the laboratory Z-axis (refer 
Fig. [9]). In Scheme 6, we have restricted the z-component 
of single-ion anisotropy to the laboratory X-Y plane. We 
have studied the variation in the molecular anisotropies 
as a function of orientation of local anisotropy in these 
schemes for ground and the excited eigenstates (Fig. 111. 
We show in Table [TTJ the Dm values for the ground and 
the excited spin states of the molecule for schemes 4, 5 
and 6. We note that when the local anisotropies are sys- 
tematically varied, there is a very large variation in the 
molecular anisotropy as a function of the local orientation 
(Fig. " 
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similar to Mn^Ac, in all the eigenstates that 
we have studied. We also note that given the orientations 
of local anisotropies, the actual molecular anisotropy val- 
ues are not very different in different spin eigenstates, 
since the energy gaps between the ground and the ex- 
cited states are small and since the spin correlations in 
these states are not significantly different. The orienta- 
tion of the local anisotropy centers for which we get the 
best agreement with experiments (Dm = -0.28 K) cor- 
responds to 6 ~ 99° [221 123]. As with Mn 12 Ac, we find 
that the laboratory frame we have chosen is very close 
to molecular axis in all the cases. In case of Fe$ cluster, 
the Z?2 symmetry commutes with the Hamiltonian in Eq. 
[2] and allows for a non-zero Em term. The variation of 
Em as a function of 9 is shown in Fig. [12] the value of 
Em for which Dm has the best fit is 0.017 K compared 
to the experimental estimate of 0.046 K obtained from 
High-frequency EPR measurements [19]. 




FIG. 9: Schematic diagram showing the directions of lo- 
cal anisotropy in Fes. The ^-component of the single-ion 
anisotropies of all the Fe(III) ions are inclined at an angle 9 
to the laboratory Z axis (Scheme 5). 




FIG. 10: Schematic diagram showing the directions of lo- 
cal anisotropy in Fes- The z-component of the single-ion 
anisotropies of all the Fe(III) ions are directed along the 
plane perpendicular to the laboratory Z axis (Scheme 6). The 
arrows indicate the Fe-O bonds on which the chosen local x- 
axis has maximum projection. 



TABLE II: Dm values of ground and excited states of Fes 
under schemes 4, 5 and 6 in K. For scheme 5, we have pre- 
sented the Dm values only for 0=98.78° for which the Dm 
value of the ground state matches with the experimentally 
observed value. 



IV. CONCLUSIONS 

In this paper we presented a general method to calcu- 
late the molecular magnetic anisotropy parameters, Dm 
and Em for single molecule magnets in a chosen eigen- 
state of the exchange Hamiltonian. Since, anisotropy is 
generally weak in SMMs compared to exchange interac- 
tion, we treat the anisotropy Hamiltonian as a pertur- 
bation over the exchange Hamiltonian. Calculation of 



State 


Dm (K) 


Scheme 4 


Scheme 5 


Scheme 6 


Ground state (S=10) 


0.6030 


-0.2892 


-0.3034 


First excited state 
(S=9) £ 3 =13.56 K 


0.5821 


-0.2783 


-0.2923 


Second excited state 
(S=9) £ 9 =27.28 K 


0.5877 


-0.2797 


-0.2952 


Third excited state 
(S=8) £ 9 =28.33 K 


0.5503 


-0.2694 


-0.2758 



7 



0.5 



0.3 



0.1 



-0.1 



-0.3 



-0.5 





G— 


-OS=10 




E— 


-0S=9 




A— 


— &S=9 




0- 


— OS=80^ 









30 



60 



90 

9 (degree) 



120 



150 



180 



FIG. 11: Variation of Dm in Feg cluster as a function of 9, 
the angle the ^-component of local anisotropy of Fe(III) ions 
makes with the laboratory Z-axis for scheme 5. 
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of molecular magnetic anisotropy compared to the ex- 
perimental values. Therefore, we focus on the molecular 
anisotropy from the single-ion anisotropics of the individ- 
ual transition metal centers in the SMM. The single-ion 
anisotropy has relativistic origin, (spin-orbit interactions 
generally dominate over dipolar interactions between un- 
paired electrons in case of transition metal ions) and are 
hence short ranged with inverse cube dependence on dis- 
tance. Therefore, we neglect interaction between spin 
moment on one ion with the orbital moment on another. 
This approximation simplifies the perturbation Hamilto- 
nian. The molecular anisotropics are computed from the 
single-ion anisotropics, using first order perturbation the- 
ory for different spin states of the SMMs. We have com- 
puted the molecular magnetic anisotropy parameters of 
Mn^Ac and Fes SMMs in various eigenstates of differ- 
ent total spin. We also studied the variation of molecular 
anisotropy by rotating the local anisotropy of the metal 
ions. In case of Mn^Ac, we find that the molecular 
anisotropy changes drastically with the local anisotropy 
direction. The Dm value is different in ground and ex- 
cited states we have computed, owing to large difference 
in spin-spin correlation values. The molecular anisotropy 
of MrinAc does not change significantly with the orien- 
tation of the local anisotropy of the core Mn(IV) ions. In 
the case of Fes cluster also, we find that the molecular 
anisotropy parameters depend strongly on the orienta- 
tion of the local anisotropy. Dm value is not very differ- 
ent in ground and excited states probably due to small 
energy gaps which implies similar spin-spin correlations. 
In case of Mn^Ac, the first order rhombic anisotropy 
term is zero due to the Did symmetry of the molecule 
while it is non-zero in Fes. The second order rhombic 
anisotropy terms commute with the molecular symmetry 
of the Mn\iAc cluster and cause tunneling between the 
states on cither side of the double potential well. Our 
method can also be extended to the calculation of these 
anisotropy constants. 



FIG. 12: Variation of Em in Fes cluster as a function of 6, 
the angle the z-component of local anisotropy of Fe(III) ions 
makes with the laboratory Z-axis for scheme 5. 

Dm and Em values assuming only dipolar interactions 
between the transition metal ions give negligible values 
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